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Abstract 

It is performed for the first time a next-to-next-to-leading order 
analysis of deep inelastic structure functions F 2 and Fl using the 
recently determined first moments of the non-singlet anomalous 
dimensions and the corresponding Wilson coefficients. From a 
comparison to BCDMS data we found a slight decrease in the 
energy scale of perturbative QCD with respect to the next-to- 
leading order analysis, which brings the scale closer to the value 
predicted from a scheme invariant approach. 



The determination of the QCD coupling constant a s at some energy is 
one of the present problems of both experimental and theoretical high energy 
physics. Deep inelastic lepton-nucleon scattering (DIS) has been one of the 
most important processes for qualitative tests of perturbative QCD and in 
particular for the determination of the energy scale A involved in a s . In the 
last few years quite exact experimental data for DIS structure functions has 
been obtained (see [jlj, for review) over a wide range of the Bjorken variable 
x. Theoretical knowledge had been restricted to next-to-leading order (NLO) 
approximation (see and its references) except sum rules |§ and the lower 
(Mellin) moment of the longitudinal structure function Fl [|J which have 
been known with next-to-next-to-leading (NNLO) accuracy. 

The purpose of this letter is to present a NNLO analysis of structure func- 
tions. The two new ingredients needed in the calculation, i.e. the anomalous 
dimensions at three loops and the Wilson coefficients of the quark and gluon 
operators at two and three loops for the transverse and longitudinal parts, re- 
spectively, were recently calculated by Larin, van Ritbergen and Vermaseren: 
The first ten 0(a 2 ) moments of the F2 coefficients 1 , and the n = 2, 4, 
6 and 8 moments for the non-singlet (NS) anomalous dimension at NNLO, 
and also the Fl coefficients at O(a^) [0]. 

Some attempts have also been made in order to include the effect of higher 
order corrections in structure functions using the information of the NLO and 
the invariance on the renormalization scheme of the physical results. Here we 
extend one of these scheme invariant approaches with the inclusion of NNLO 
terms. 

At three loops the Mellin moments of structure functions at two different 
scales Q 2 and Ql are related byQ 
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1 The corresponding kernels into Mellin convolutions with parton distributions were 
calculated by van Neerven and Zijlstra |^| 

2 Through this work we use the short notation a — a s (Q 2 ) / 'Air and «o = o; s (Qq)/4:Tt, 
and the standard definition for moments M^ n = dxx n ~ 2 Fi(x) . For more details on the 
notation see ref. [01 and references therein 
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where j extends to singlet and non-singlet contributions. 

The functions d], Q$ j and are related to the one- loop, two- loop [|9f 
and three-loop ]7] anomalous dimensions ('yW™, % = 0, 1, 2). In the non-singlet 
case one has 
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The quantities B££ (S^) and B)*» (B£») are the NLO (NNLO) Wilson 
coefficients of F 2 and Fl respectively. References for the two-loop coefficients 
can be found for example in [0, while the new three- loop results are in || 
and 0. In eq. (|) R { ^ S = Bf£ s /B^ 8 for i = 2, 3. fa are the coefficients 
of the renormalization group beta function. For the scheme dependent fa we 
use the MS 1 expression given in ref. |§ 

In eq. (|TJ) the moments M 2 (Ql) and M L (<5o) are related at NLO by 

MIM) = Ml n (Ql)a B^(l + a [R < Si-B^} 

+ a 0[ K L,n ~ B 2,n ~ &2,n \ K L,n ~ B 2,n )\ 

and the three-loop coupling constant is obtained by solving the transcenden 
tal equation 
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The only piece unknown in the perturbative calculation is the moment of 
F 2 at some energy scale Qq. In the non-singlet case it is parametrized by 
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where A, a, (3 and 7, in addition to the energy scale A, have to be obtained 
from fits to experimental data. 

The expressions for the moments M 2 n (Q 2 ) and M^^Q 2 ) in the scheme- 
invariant (SI) approach can be found in 0, where the equations for the SI 
coupling constants ajj.„, for i=l (NLO), were also given. At three loops they 
obey the following equations: 



A) In 
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In the two-loop approximation the couplings a^J obey similar equations to 

eq. (|j) but without the term linear in a^J . 

With the statistic of present experiments it has been shown that structure 
functions (F) can be reliably reconstructed from their moments with the help 
of orthogonal Jacobi polynomials Q% (x) (see |H| and references therein) 

■^V max 

F(x,Q 2 ) ~ x a {l-xfY. a k (Q 2 )Qf(x), (9) 

k=0 
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where the coefficients are connected to the moments of F through 



MQ 2 ) = £cf (a,/?)M j+2 (Q 2 ). (10) 

3=0 

This method has been used as a simple and fast alternative (in computer 
consuming time) to solve the evolution equations, but here it is specially 
useful because only the moments are known instead of the explicit analytical 
result in x. 

Notice that for the non-singlet anomalous dimensions and Fl Wilson 
coefficients only the first four even moments have been calculated (see above) 
while the polynomial reconstruction requires also the odd ones and at least 



8 in total in order to get an accurate result |14j]. We solve this problem 
by obtaining the values for n = 3, 5 and 7 from smooth interpolation (see 
fig. 1). This interpolation is based on theoretical investigations (see |15], [16|1 ) 
where it was showed that the n-dependent functions in the Wilson coefficients 
and anomalous dimensions can be transformed to polygamma-f unctions (and 
some others associated to it) which are continuous in its arguments. Thus, 
the analytical extension from even values of the argument to odd ones (see 
[ |j~5|j ) and moreover to non integer ones (see [[RJ) can be done without any 
problem. 

The three-loops anomalous dimension can also be extrapolated to n > 8 
if one assumes that the behavior observed for the first eight moments, which 
is similar to the two-loop anomalous dimension (see fig. 2), remains at higher 
n. In that case one can use the relation 

(2)2 

7$r = 7^ (id 

Tns 

to extend the calculation to higher n (dashed line in fig. 1). 

There is not, however, an analogous relation to extrapolate the Wilson 
coefficients, of which we only know the values at n=2,...,10 for F 2 and n = 2, 

8 for Fl. Thus the sum in eq. (^) was restricted to N max =6, i.e. moments 
from n=2 to n=8. 

The inclusion of target mass (TM) corrections is not possible at NNLO, 
because one would need to know the moments for n > 8 (see eq. (23) from 
H). Also, we did not consider TM effects in the rest of the analysis in order 
to make a meaningful comparison between the results. 
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We have performed a series of non-singlet fits to F2 BCDMS deuterium 
ll|] and proton data, and also to SLAC reanalized R = Fl/2xF\ data 



from deuterium [ID|. Structure functions were calculated at LO, NLO, NNLO 
and in both scheme invariant, NLO (SI-N) and NNLO (SI-NN), approaches 
for comparisonP]. Throughout this work we have considered nf = 4 in the 
theoretical expressions, and the reference point in eq. (|5D was fixed to Q$ = 
50 GeV 2 . 

In the analysis of R we have fitted simultaneously F2 in order to fix as 
much as possible the scale A with the help of the more precise BCDMS data. 

The analysis has been restricted to data with Q 2 > 5 GeV 2 to avoid the 
kinematic region where higher twist effects could affect F2 because they are 
not under theoretical control. In R this cut excludes the region affected by 
corrections of order 1/Q A or higher. Also we analyze the region x > 0.35 
where the influence of gluons can be neglected and the non-singlet evolution 
is justified. In addition, it has been used F 2 data with y > 0.16, being y 
the fraction of energy in the laboratory system transferred to the nucleus, in 
order to exclude a region of large systematic errors. In practice it produces 
a different Q 2 cut for each x bin which excludes even more the region (high 
x and low Q 2 ) where power corrections could affect. 

Table 1 summarizes the results obtained. In the fits to only F2 we find 
that there is a small difference between \ 2 obtained from different theoreti- 
cal predictions, which shows that data does not discriminate between them. 
However, in the analysis of the more precise proton data, the LO fit gives 
the worst agreement. 

There are also differences in the central value of A of the fits. The scale 
obtained in the SI-N analysis, even though within the errors, is significantly 
smaller than the NLO result (12 % for deuterium and 9 % for proton). This 
decreasing was already observed in earlier similar analyses (see for example 
ref. @ and references therein) and interpreted as the result of the influence 
of higher order corrections which are better taken into account by the scheme 
invariant treatment^. Now, one can test this supposition by comparing the 

3 The results at LO, NLO and SI-N for deuterium here are different from those in a 
similar analysis in reference Q due to the inclusion of target mass corrections in . We 
have also checked that the using of a higher N max in eq. (||), when it is possible, does not 
affect significantly the results 

4 We do not clarify the concrete SI procedure here. This may be Stevenson's PMS 
]22j, Grunberg's EC method [E3 and other procedures (see for example [0, E4j and its 
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results from the analysis NNLO and SI-N. One see in table 1 that the value 
of A at NNLO is lower than at NLO even though the decreasing is not so 
pronounced as it is in the SI-N case. We think however that it shows the 
trend predicted by SI-N. This behavior has also been observed in another 
independent analysis pO . 



We have also fitted F 2 in the SI-NN case. In this case equations @ 
above and also eq. (16) and (17) in |2j have been used in the calculation. 
The results presented in table 1 show that there is no appreciable variation 
in A with respect to the standard NNLO approximation. 

For the combined fits to F 2 and R (without twist corrections) the situation 
is unfortunately more complicated. Some of the predictions are fitting better 
R than F 2 , giving higher A values that can not be compared. Notice for 
example the poor partial fit to F 2 in the SI-N analysis. However, if target 
mass corrections are included in R, one gets a very good fit to F 2 (see for 
example |2|). It shows that these corrections, which are very large in R 
[@, ^l|] , are essential for a meaningful combined analysis. 

Comparing the results obtained from the NLO and NNLO analyses, one 
can also notice, as in the fits to only F 2 , a small decreasing in A in the NNLO 
case. 

In the SI-NN fits we were not able to converge to a reasonable value of 
X 2 . This could be due to large a 3 corrections in Fl, or again an effect due to 
the absence of TM corrections. Also it is possible that the use of N max = 6 
to reconstruct Fl from the moments is not sufficient in this case. This point 
will be clarify when a 3 correction to Fl are known at larger values of n. 

Finally, trying to improve the agreement with R data we have included 
twist-4 corrections in the calculation of Fl. That contribution has been 
simply parametrized by |17j 

M ^-4 (g2) = 8 ^ M2n(g2)) (12) 

where the energy scale k is the only unknown parameter which has to be 
determined from data. 

In all cases we found a significant improvement of the agreement with 
R data. The scale k is higher than in previous analysis because in the 
present work it contains part of the effect due to target mass corrections. 



references) , which seems to lead to similar conclusions (sec [nM [L9| ) 
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The value of A from the NNLO analysis is lower than the NLO result, as in 
the others cases analyzed above, and equal to the SI-NN result. For the scale 
k, this behavior is obtained only in the proton case. 

In conclusion, it has been performed for the first time a NNLO analysis of 
F 2 and R data. We found a slight decrease of A when the NNLO expression 
for the moments is used, which brings the result closer to scheme invariant 
predictions. This behavior is also in agreement with similar analyses of R- 
ratio in e + e~ annihilation, the r-lepton decay rate (see [[L8[]) and the Gross- 
Llewellyn Smith sum rule in |Il9| . 
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Figure captions 



1. Three-loop coefficient of the non-singlet anomalous dimension in the 
MS renormalization scheme. 

2. One, two and three-loop coefficients of the non- singlet anomalous di- 
mension in the MS renormalization scheme normalized to the value of 
the first moment n = 2. 



Table captions 



1. Results of non-singlet fits to F 2 proton [11]] and deuterium [T2j data 
from BCDMS, and to R data from SLAC [|l(|. Errors given are statis- 
tical. 
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Deuterium 




(MeV) 


k (MeV) 


X 2 (F 2 )/dof 


X 2 (^)/dof 


F 2 


LO 


182 ±32 




63.5/65 






NLO 


182 ±30 




62.7/65 






SI(N) 


159 ±25 




62.4/65 






NNLO 


168 ±27 




62.5/65 






SI(NN) 


164 ± 26 




62.4/65 




F 2 and R 


LO 


223 ± 35 




65.1/65 


91.5/43 




NLO 


235 ± 34 




65.5/65 


90.9/43 




SI(N) 


238 ± 28 




70.8/65 


66.3/43 




NNLO 


218 ±27 




65.5/65 


79.9/43 


F 2 and R 


LO 


181 ±31 


207 ± 13 


63.5/65 


41.4/43 


(Including twist-4) 


NLO 


180 ±29 


198 ± 14 


62.7/65 


41.5/43 




SI(N) 


157 ±25 


184 ± 17 


62.4/65 


41.4/43 




NNLO 


159 ±21 


200 ± 12 


62.8/65 


44.6/43 


Proton 












F 2 


LO 


171 ±27 




51.2/60 






NLO 


175 ±26 




47.6/60 






SI(N) 


159 ±23 




46.0/60 






NNLO 


168 ±25 




46.4/60 






SI(NN) 


169 ±25 




45.8/60 




F 2 and R 


LO 


200 ± 30 




52.8/60 


92.3/43 




NLO 


222 ± 29 




50.0/60 


82.8/43 




SI(N) 


231 ±25 




54.5/60 


59.2/43 




NNLO 


209 ± 23 




48.8/60 


79.7/43 


F 2 and R 


LO 


168 ±27 


205 ± 14 


52.0/60 


43.0/43 


(Including twist-4) 


NLO 


175 ±26 


195 ± 14 


48.2/60 


43.5/43 




SI(N) 


158 ±23 


179 ± 17 


46.6/60 


42.7/43 




NNLO 


158 ±22 


189 ± 14 


46.5/60 


47.8/43 



11 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405290vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405290vl 



